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John  Kitchin,  Naftali  A.  Langberg,  and  Frank  Proschan 

ABSTRACT 

We  construct  a  new  estimator  for  a  continuous  life  distribution  from  incomplete 
data,  the  Piecewise  Exponential  Estimator  (PEXE) .  We  show  that  the  PEXE  is  strongly 
consistent  under  a  mild  restriction  on  the  distribution  of  the  censoring  random 
variables  (possibly  non- identical  and  non - cont inuous ) .  Then  we  consider  the  Product 
Limit  Estimator  (PLE) ,  introduced  by  Kaplan  and  Meier  (1958) .  We  prove  the  strong 
consistency  of  the  PLE  under  a  mild  regularity  condition  on  the  distributions  of  the 
censoring  random  variables.  This  result  extends  previous  ones  obtained  by  various 
researchers.  Finally  we  compare  the  new  PEXE  and  traditional  PLE. 
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1.  Introduction  and  Summary. 


Let  Xp  X^,  ...  be  independent  identically  distributed  (i.i.d.)  lifelengths  with 
a  cannon  continuous  life  distribution  G  and  let  L^>  •• •  be  nonnegative  censoring 

random  variables  (r.v.'s)  not  necessarily  continuous  or  i.i.d.  We  assume  that 
[Xj,  L^J,  CX2,  L2D,  ...  is  a  sequence  of  independent  pairs  defined  on  a  comnon 
probability  space  (n,  B,  P ) .  Without  loss  of  generality  we  assume  that  for  each 
in  c  fl  the  sequence  X^w),  X2(u>)>  ...  consists  of  distinct  positive  real  numbers. 

Let  I  denote  the  indicator  function,  let  T  ~  min{X  ,  L }  be  the  removal  time  of 

4  4  4 

the  qth  item  on  test,  and  let  r  =  I(X s  L  )  denote  the  cause  of  removal  of  the  qth 

4  4  4 

item,  q  *  1,  2,  ...  .  We  consider  the  problem  of  estimating  G  =  1  -  G,  the  under¬ 
lying  survival  function,  from  CT^,  ,  ...»  [Tr,  cn-i»  n  *  1,  2,  ...  .  This  is  the 
incomplete  data  problem  as  formulated  by  Kaplan  and  Meier  (1958) .  We  note  that  our 
assumptions  on  the  lifelengths  and  censoring  r.v.’s  are  less  restrictive  than  those 
assumed  in  the  theory  of  competing  risks,  the  theory  of  life  tables,  and  the  usual 
treatments  in  biostatistics.  Traditionally  one  or  more  of  the  following  assumptions 
have  been  made:  i.i.d.  censoring  r.v.’s,  continuous  lifelengths,  continuous  censoring 
r.v.'s,  and  the  independence  of  the  lifelength  and  the  corresponding  censoring  r.v.  of 
each  item  on  test.  [Kaplan  and  Meier  (1958),  Breslow  and  Crowley  (1974),  Peterson 
(1977),  Langberg,  Proschan,  and  $iinzi  (1980) J. 

In  Section  2  we  construct  a  new  estimator  for  the  underlying  survival  function  G, 
the  Piecewise  Exponential  Estimator  (PEXE) ,  and  denote  it  by  En(t) .  In  Section  3  we 
use  a  theorem  proven  in  Section  4  to  show  that  the  PEXE  is  a  strongly  consistent 
estimator  of  IT  under  mild  regularity  conditions.  In  particular,  we  obtain  the  strong 
consistency  of  the  PEXE  when  Xq,  Lq  are  independent  r.v.'s,  q  ■  1,  2,  ...  and  under 
most  of  the  "traditional"  assumptions  discussed  in  the  previous  paragraph.  In  Section 
5  we  consider  the  Product  Limit  Estimator  (PLE) ,  introduced  by  Kaplan  and  Meier  (1958) 


ifljjllrtSnr  iUM 


v  r^_ i,' :  j 
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that  serves  as  the  principal  nonparametric  estimator  to  date  of  the  survival  function 
G.  We  show  that  under  a  variety  of  conditions  the  PLE  is  a  strongly  consistent 
estimator  of  S’.  These  results  extend  those  obtained  by  Peterson  (1977) ,  and  Langberg, 
Proschan,  and  Quinzi  (1980) .  Finally,  in  Section  6  we  compare  the  new  PEXE  and  the 
traditional  PLE. 
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2.  Piecewise  Exponential  Estimator. 

In  this  section  we  introduce  a  new  estimator  for  a  continuous  life  distribution 

from  incomplete  data:  the  Piecewise  Exponential  Estimator  (PEXE) . 

We  approach  the  incomplete  data  problem  from  the  viewpoint  of  reliability  and 

life  testing.  An  item  at  age  zero  is  placed  on  life  test.  It  eventually  leaves  the 

test  either  because  it  fails,  yielding  a  complete  life  length,  or  because  it  is 

withdrawn  while  still  functioning,  yielding  a  censored  lifelength.  Thus  if  an  item 

on  test  is  an  observed  failure,  5^  *  1,  and  if  it  is  a  withdrawal,  *  0, 

q  »  1,  ...,  n.  Starting  with  a  sample  of  initial  size  n,  the  number  of  items  at 

n  n 

time  t  is  denoted  by  N  (t)  =  £  I  (T  >  t) .  Let  x(n)  =  £  denote  the  number  of 

n  q*l  "  q*l  H 

observed  failures  and  let  2^  <  ...  <  Zn;T(n)  denote  the  consecutive  observed 

failure  times,  with  Zn>g  =  0. 

Let  R(t)  *  -in  S"(t)  denote  the  hazard  function  of  the  life  distribution  G, 
t  e  CO,  Sup{u:  C(u)  >  0».  On  the  interval  ( Zn;q_1^  Zn;qj,  we  estimate 

CZ»;,  -  Zn:q-irl«W  '  RBn:q-X))  *  V,  *  C  f"  V"1*"'1'  ^  ^  °f 

n:q-l 

observed  failures  per  unit  time  in  the  interval  (Zn.q_i»  zn:q^’  q  * •  * 5  * 

Z 

We  note  that  Jn:q  N  (u)du  is  the  total  time  on  test  in  the  interval 

z  “ 

^n:q-l 

(Z  ,,  Z  ],  q  »  1,  ...»  T(n).  These  hazard  slope  estimators: 
v  n:q-l*  n:q  ’  n  ’ 

r  ,  q  *  1,  ...,  x(n),  define  a  piecewise  linear  estimator  of  the  hazard  function  R, 
n,q 

which  in  turn  leads  to  a  piecewise  exponential  estimator  of  the  underlying 
survival  function  C,  given  explicitly  by  the  following  definition. 

Definition  2.1.  For  x(n,  u»)  i  1  let  A„  „(w)  * 

""  ~  '  n>H 

(Z„.,»  -  Z  _  !  (o»))r„.n(u>),  q  -  1,  ...,  x(n,  <o),  u  e  n.  Then  the  Piecewise 
n.q  n*q“X  n.q 

Exponential  Estimator  (PEXE)  of  the  survival  function  C,  denoted  by  E^(t,  u),  is 

?s*ual  V*  1  on  the  set  {x(n,  u>)  *  0  or  t  e  (-•,  0J),  is  equal  to 

exp {-qiLn  jCw)  -(t  -  Zn.q.jUDJx^qCw)}  on  the  set  (x(n,  w)  i  1,  t  6 
3*1  *  x(n,w) 

Zn.q(w)),  q  -  1,  ...,  t(n,  o»)>,  and  is  equal  to  exp  -  £  *n>j>  on  the  set 
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For  the  sake  of 


1  (Zn:TCn,»)W>  *>  ' 


simplicity  we  supress  the  argument  w  in  u) 


In  this  section  we  use  Theorem  3.2,  stated  below  and  proven  in  Section  4,  to 
obtain  the  strong  consistency  of  the  PEXE  under  various  conditions.  We  need  a 
definition  and  some  notation. 

Definition  3.1.  Let  K  be  a  function  defined  on  (-«,  •).  We  say  that  K  is  a 
subdistribution  function  (s.d.f.)  if  K  is  nondecreasing,  right  continuous,  and  assumes 
values  in  CO,  lj. 

For  a  s.d.f.  K,  let  K(x)  *  lim  K(y)  -  K(x)  be  the  subsurvival  function 

y-*» 

corresponding  to  K,  let  C(K)  be  the  set  of  all  continuity  points  of  K,  and  let  a(K)  = 
sup{t:  K(t)  >  0>. 

We  are  now  ready  to  state  Theorem  3.2. 

Theorem  3.2.  Assume  the  following: 

iH 

(3.1)  There  is  a  s.d.f.  F(t)  such  that  li®  n  £  I (T  >  t)  «  F(t)  for  t  €  C(F), 

■  q»l  ^ 

and 

(3.2)  There  is  a  s.d.f.  F(t,  1)  such  that 

-1  ^ 

ii':  ~  £  P{T  St,  5  ■  1}  ■  F(t,  1)  for  t  e  (-»,  •). 
i  r*»  q»l  4  4 

Then  there  is  a  set  ^  e  B,  PCn^  «  1,  such  that  for  all  u  £ 

t  r 

lim  Fn(t)  *  exp{-/cF(u)J  dF(u,  1))  for  t  e  CO,  «(F(*,  1))). 
n+®  0 

Now  we  use  Theorem  3.2  to  obtain  the  strong  consistency  of  the  PEXE  under 
various  conditions.  For  simplicity,  we  denote  throughout  a(F(*,  1))  by  a,. 

First  we  prove  the  strong  consistency  when  the  lifelength  of  each  item  on 
test  and  its  censoring  r.v.  are  independent. 


Theorem  3.3.  Assume  the  following: 

(i)  The  r.v.'s  X ,  L  are  independent  for  q  *  1,  2, 

4  M 

and 

l  n 

(ii)  There  is  a  s.d.f.  H(t)  such  that  lim  n  x  l  P{L  a)*  H(t)  for  t  e  C(H). 

n*»  q=l  ** 

Then  there  is  a  set  P(n^}  »  1,  such  that  for  all  u  e 

lim  ^(t)  *  C(t)  for  t  e  CO,  a.). 
n*»  1 

t 

Proof.  Let  F(t)  -  C(t)IT(t)  and  F(t,  1)  =  /IT(u)dG(u)  for  t  e  CO,  •).  To 

0 

obtain  the  desired  result  it  suffices,  by  Theorem  3.2,  to  verify  Conditions  (3.1), 

(3.2)  and  to  show  that: 

t  7 

(3.3)  tT(t)  *  exp{-/CF(u)J  dF(u,  1)}  for  t  e  CO,  a,). 

0  1 


First  we  verify  Conditions  (3.1),  (3.2).  Let  t  e  (-«,  •)  and  n  »  1,  2,  ... 
Then  by  Assumption  (i) : 

.n  .  n 

l 

q= 


(3.4)  V»umll  P(T  >  t>  -  (TWCn-1  £  P(L  >  t»,and 
r**«  q-1  "  q=i  1 


-ln  t  1  n 

(3.5)  li~n  l  P  T  *  t,  $  -  1}  -  /Cn‘x  £  P{L  >  u}JdG(u). 
n**»  q»l  4  4  0  4 

By  Assumption  (i)  and  the  definition  of  F,  C(F)  »  C(H).  Thus  Condition  (3.1) 
follows  by  Assumptions  (i),  (ii),  and  Equation  (3.4). 

Condition  (3.2)  follows  from  Assumptions  (i),  (ii),  Equation  (3.5),  and  the 
dominated  convergence  theorem. 

Finally  we  verify  Equation  (3.3).  By  the  definitions  of  F(t)  and  F(t,  1), 
t  ,  t  , 

^CF(u)]  dF(u,  1)  *  Jc£(u)]  dG(u)  for  t  c  CO,  a(F)).  Consequently  Equation  (3.3) 
follows  by  the  continuity  of  G .  | | 
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Next  we  obtain  from  Theorem  3.3  a  corollary  of  practical  importance. 

Assume  that  most  of  the  "traditional"  conditions  presented  in  Section  1  hold 
(but  not  necessarily  the  continuity  of  the  censoring  r.v.'s).  Then  Assumptions  (i), 
(ii)  of  Theorem  3.3  hold  and  we  obtain: 

Corollary  3.4.  Assume  that  the  r.v.'s  X  ,  L  are  independent,  q  =  1,  2,  . .., 
and  that  the  r.v.'s  L^,  q»l,  2,  ...»  are  i.i.d.  Then  there  is  a  set  e  B,  P{«1}  =  1 
such  that  for  all  u  c 

lim  F  (t)  *  C(t)  for  t  e  [0,  o.) . 

Now  we  derive  the  strong  consistency  of  E"n(t)  when  the  censoring  r.v. 's  are  i.i.d 

but  not  necessarily  independent  of  the  corresponding  lifelengths. 

Theorem  3.5.  Assume  that  the  pairs  LX^,  L^],  q  »  1,  2,  ...  are  i.i.d.  Then 

there  is  a  set  fl.  e  B,  P{n.)  ■  1,  such  that  for  all  wen.,  lim  EL Ct)  *  (T(t)  for 

1  n^» 

t  €  [0,  Oj)  iff 

(3.6)  P{LX  >  t|Xx  -  t)  *  P^  >  t|Xx  >  t}  for  t  «  [0,  . 

Proof.  To  obtain  the  desired  result  it  suffices  to  show  that  (a)  Conditions 
(3.1),  (3.2)  hold  and  that  (b)  Conditions  (3.3)  and  (3.6)  are  equivalent. 

(a)  Let  F(t)  *  PO^  s  t)  and  F(t,  1)  *  P{Tj  s  t,  ^  »  H  t  «  (-*»,  ») .  Then 
Conditions  (3.1),  (3.2)  hold  trivially. 

(b)  Let  A  be  a  Borel  set  contained  in  CO,  a^).  Then: 

(3.7)  P{TX  e  A,  -  1}  -  /  P^  >  u|Xx  *  u}dG(u). 

A 

First  we  show  that  Condition  (3.3)  implies  Condition  (3.6).  By  the 
continuity  of  G  and  by  Equation  (3.7): 
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t  ,  t  i 

/[{T(u)  J’aG(u)  *  -tnG(t)  *  /[F(u)  J' AdF(u,  1) 

0  0 

1  1 

-  /[F(u)]P{L,  >  u|X.  =  u)dG(u)  for  t  c  [0,  a.). 

0  x  i  i 

Consequently  Condition  (3.6)  follows  since  5  a(G). 

Finally  we  show  that  Condition  (3.6)  implies  Condition  (3.3).  By  Condition 
(3.6)  and  Equation  (3.7),  we  have: 

t  1 

(3.8)  E(t,  1)  *  J [ ff(u) 3  "T (u)dG (u)  for  t  e  [0,  a,),  and 

0  1 

t  ,  t  , 

(3.9)  J[F(u)j  dF(u,  1)  =  /(^(u) j'aG(u)  for  t  e  [0,  a.). 

0  0  1 

Consequently  Condition  (3.3)  follows  from  Equation  (3.9)  and  the  continuity  of  G.  | | 
Note  that  E"n(t)  is  a  strongly  consistent  estimator  of  the  survival  function  (T 
whenever  Conditions  (3.1)  through  (3.3)  hold.  Finally  we  provide  an  example  where 
Conditions  (3.1)  through  (3.3)  hold.  First  we  need  the  following  definition. 

Definition  3.6.  [Marshall  and  Olkin  (1967) J.  Let  X^,  X2,  and  X^  2 
nonnegative  real  numbers,  x^  +  X2  +  x^  2  >  0.  Then  the  random  pair  CU,  V]  with 
nonnegative  components  has  a  Marshall  Qlkin  Bivariate  Exponential  Distribution 
(MOBVED)  with  parameter  [X^,  x2,  x^  2J  if  for  all  t,  s  e  [0,  «): 

P{U  >  t,  V  >  s)  *  exp{-A^t  -  x2s  -  X^  2  max(t>  s)}. 

We  now  construct  the  example. 

Example  3.7.  Let  X^,  x^  2,  y^,  y2,  ...  be  nonnegative  real  numbers, 

A  ■  X,  +  X,  ,  ♦  y  >  0,  q  *  1,  2,  . . .,  and  let  [X  ,  L  j,  q  *  1,  2,  ...,  be  a  sequence 
q  l  q  q  q 

of  random  pairs  having  MOBVED’ s  with  parameters  CX^,  Yq>  Aj  2J  respectively.  Assume 

that  lim  *  y  «  (0,  »)>  and  that  X.  +  X.  ,  +  y  >  0.  We  show  that  under  these 
n*a,  n  l 

assumptions,  Fn(t)  is  a  strongly  consistent  estimator. 
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Note  that  for  t  c  [0,  »)  awl  q  *  1,  2, 

P{Xq  >  t>  *  expMAj  +  x1>2)t),  PtTq  >  t>*  exp{-Aqt},  and 

WTq  >  t,  tq  -  1)  •  C»J  *  »l,2UqlexP{-*qt>' 

Let  F(t)  -  exp(-At)  andF(t,  1)  *  (xx  +  x1>2)A'1exp{-At},  t  e  CO,  «).  Clearly 
Conditions  (3.1)  through  (3.3)  hold.  Thus,  by  Theorem  3.2,  En(t)  is  a  strongly 
consistent  estimator  of  G  for  t  e  CO,  ») . 

In  particular,  if  the  random  pairs  [Xq,  LqJ,  q  *  1,  2,  ...,  are  i.i.d.  with  a 
MOBVED,  then  E^(t)  is  a  strongly  consistent  estimator  of  S’  for  t  e  [0,  ®) . 
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4.  Proof  of  Theorem  3.2. 


In  this  section  we  present  a  proof  of  Theorem  3.2.  Let 
.  n 


F_(t,  o>)  -  nll  ICT>)  >  t)  and  F(t,  1,  «) 
n  q*l  **  n 


n"  £  I(T  W  ^  t,  £ (<d)  *  1),  n  =  1,  2,  ...»  t  c  (-«,  •),  u  e  fl. 
q»l  4  4 

For  the  sake  of  simplicity  in  notation  we  supress  throughout  the  argument  u. 

t  , 

Note  that  the  continuity  of  the  function  exp(-/CF(u) J  dF(u,  1)}  in  CO,  o1) 

0  1 

follows  from  the  continuity  of  F(t,  1),  which  in  turn,  follows  from  the  continuity 
of  G.  Thus  to  prove  the  result  of  Theorem  3.1,  it  suffices  to  show  by  a  standard 
argument  [Chung  (1974),  pp.  132-133]  that: 


(4.1)  lim  Eft)  ■  exp{-J[F(u)]'AdF(u,  1)),  w.p.l.  for  t  e  [0,  a,). 

n-*“  n  0  1 

First  we  prove  that  to  obtain  Statement  (4.1),  it  suffices  to  show  that: 
t  t  t  , 

(4.2)  lim  J[F(u)]‘AdF_(u,  1)  *  /[F(U)3  dF(u,  1),  w.p.l.  for  t  «  [0,  o(F)) 

m- 0  n  n  0 

Then  we  complete  the  proof  of  Theorem  3.2  by  verifying  Statement  (4.2). 

We  need  some  more  notation  and  two  lemmas.  Let  a(n,  t)  ■ 
maxfk:  k  *  0,  ...»  x(n),  Zn<k  <  t),  t  e  [0,  aj),  n  «  1,  2,  ....  and  let  CxJ 
denote  the  largest  integer  less  than  or  equal  to  x,  x  e  (-“,  •) . 


Lenina  4.1.  Assume  that  Conditions  (3.1),  (3.2)  hold.  Then 


(4.3) 


(4.4) 


lim  Fn(t)  -  F(t),  w.p.l.  for  t  e  C(F), 
n**> 


lim  Fn(t,  1)  *  F(t,  1),  w.p.l.  for  all  t  e  [0,  •). 
n*» 
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Proof.  To  prove  Statements  (4.3)  and  (4.4)  it  suffices  to  show  by  Conditions 
(3.1),  (3.2)  that  for  all  t  e  (-•,  «): 

-l  n 

limEF  (t)  -  n  J  P{T  >  t}j  *  0,  w.  p.  1.,  and  that 
n*~  q-1  q 

-1  n 

limCF  (t,  1)  -  nx  l  F{T  s  t,  5  -  1)3  ■  0,  w.p.l. 
rH°°  q»l  q  q 

The  preceding  two  statements  follow  by  the  strong  law  of  large  numbers .  | | 
Lemma  4.2.  Assume  that  Conditions  (3.1),  (3.2)  hold.  Then 

(4.5)  li5Lzn:T(n)  >  w.p.l.  for  t  c  CO,  c^) . 

Proof.  Let  lim  Z _ =  a.  For  a  *  »,  Statement  (4.5)  follows  trivially. 

n*  '■n-' 

Assume  a  <  »;  then  by  Statement  (4 .4) : 

*  .1  .1  m 
F(s,  1)  *  lim  m  F  (s,  1)  s  limm  1  J  p{E  -  1} 

q»i  q 


s  F(a,  1)  for  all  s  c  CO,  «•),  where  (m)  is  a  subsequence  of  the  positive 


integers . 


Thus  a  *  a^.  Consequently  Statement  (4.5)  follows.  || 

We  show  now  that  to  prove  Statement  (4 .1) ,  it  suffices  to  verify  Statement 
(4.2). 

Lemma  4.3.  Assume  that  Conditions  (3.1),  (3.2)  hold.  Then: 

t  , 

(4.6)  lim |  in  IT  (t)  ♦  J[F  (u)3"AdF  -  0,  w.p.l.  for  t  e  CO,  a.). 

|>*»  0  n>*  A 

Proof.  Let  fl 2  «  anc^  let  1  «  CO,  “i ) •  Then  by  Statements  (4.3) 

through  (4.5)  for  every  u  e  n,,  there  is  a  positive  integer  n(«)  such  that 

t  Zi 

Fn(t)  >  0,  t(n)  a  1,  ^(Fn(u)J  dFn(u,  1)  <  »,  and  Zn;t^  >  t  for  n  2  n(«).  Let 
u  e  ft2.  Then: 


3  2 


(4.7)  [nfn(Zn.q.l)j-1  4  (Zn;q  -  ^n!q.1)rn(q  4 

[nF_(Z  )]_1  for  q  *  1,  t(n),  n  a  n(w). 

n  n.q 

By  Inequality  (4.7): 

a(n,t)  .1 

-tnEn(t)  a  I  tnFn(Zn  ,j)j  i 
q*l 

t  n  .1 

J[F  (u)3  dF  (u,  1)  -  CnF  (t)j  for  n  a  n(«). 

0  "  n 

Consequently: 

t  i 

(4.8)  Iun[tnIL(t)  +  J[F  (u)J  dFn(u,  1)3  ^  0 . 

n*»  0 

Further,  by  Inequality  (4.7): 

a(n,t)+l  it  , 

-*riTn(t)  4  {  cnrnczn  J3-1  4  fCFn(u))  ^Cu,  1) 

q*l  u 

+  tnFn^n:a(n,t)+l^  . 

By  Condition  (3.2)  there  is  a  6  e  (0,  -  t)  and  a  positive  integer  nx  such 

a(n,  t)  +  1  s  nF(t  +  6,  1)  for  n  a  n^  Thus: 

tnFntZn:a(n,t).l)J'1  s  Cn(1  ’  n'1‘a(n"  tJ  *  1))J‘1 
s  [n(l  -  F(t  +  6,  DU'1  *  CnF(t  +  6,  DU*1  for  n  a  i^. 
Consequently: 

t  i 

(4.9)  liin  ClnE_(t)  ♦  /CFn(u)j  dFn(u,  1)3  a  0. 

n*»  0 

Statement  (4.6)  follows  now  from  Inequalities  (4.8)  and  (4.9).  || 


Let  0  «  Sup{F(t,  1):  t  e  [0,  «»)},  t  c  [0,  a(F)),  and  b(n,  t) 


max{-r(n),  a(n,  t)  +1}.  Assune  that  the  r.v.'s  Lp  ....  L2,  ...  are  i.i.d.  Then  by 
considering  the  cases  F(t,  1)  ■  0  and  F(t,  1)  <  0,  it  follows  that 
Um  F^CZn.^^n  >  0.  Thus  if  the  r.v.'s  L^,  L2,  ...  are  i.i.d.  Statement  (4.6) 
holds  for  t  e  CO,  o(F)). 

We  complete  the  proof  of  Theorem  3 . 2  by  verifying  Statement  (4.2). 

Lenina  4.4.  Assume  that  Conditions  (3.1),  (3.2)  are  satisfied.  Then 


Statement  (4.2)  holds. 


t  e  CO,  a(F)).  Upon  integration  by  parts; 


1),  u  e  (-»,  •),  and  let 


^  _1  t  -l 

/[F_(u)j  AdF  (u,  1)  =  -jF(u-,  l)d[F  (u)j'A 
0  n  n  0  n  n 

♦  CFn(t))'1Fn(t,  1)  -  CFn(0)j'1Fn(0,  1)  - 

*/{F  (u-,  1)  -  F(u,  l))d(F hWJ^./Vu,  l)d[F  (tt)j*1 
0  11  “  0 


+  1)  -  [Fn(0)j'1Fn(0,  1). 


Thus,  by  the  continuity  of  F(«,  1),  and  upon  integration  by  parts: 


(4.10) 


t  1  t  1 

J[F(u)JdF(u,  1)  -  /(F  (u)J~oF(u,  1)  + 
0  n  0  n 


CFn(t)r1{Fn(t,  1)  -  F(t,  1))  -  CFn(0)]'1{Fn(0,  1)  -  F(0,  1)> 


t  , 

-  ^(Fn(u-,  1)  -  F(u,  l)>d[Fn(u)j  \ 

Next  note  that  by  Statement  (4.3),  the  continuity  of  F(u,  1),  and  the  dominated 
convergence  theorem: 

t  1  t  | 

(4.11)  lim  /[F_(u))‘AdF(u,  1)  -  /(F(u))'AdF(u,  1),  w.p.l., 

Tt*~  0  0 


I  *• 

t’v.t  by  Statements  (4.3),  (4.4): 

(4.12)  limCFn(t)j'1{Fn(t,  1)  -  F(t,  1)>  »  0,  and 

limCFL(O) j"*{F  (0,  1)  -  F(0,  1)>  =  0,  w.p.l., 
n-*x>  n 

and  that  by  the  continuity  of  F(* ,  1): 

(4.13)  lira  Siq3{|Fn(u-,  1)  -  F(u,  1)|,  u  «  [0,  tj>  *  0,  w.p.l. 
n*» 

Consequently  the  desired  result  follows  by  Statements  (4.10)  through  (4.13).  || 
Finally, assume  that  the  r.v.'s  L^,  I^.  •••»  are  i.i.d.  Then  by  the  remark 
following  Lemma  4.3,  and  1/  Lemma  4.4,  Statement  (4.1)  holds  for  t  e  CO,  a(F)). 
Consequently  if  the  r.v.'s  L^,  L2,  ...  are  i.i.d.  then  Corollary  3.4,  holds  for 
t  e  [0,  a(F)) . 
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In  this  section  we  first  prove  that  as  n  +  »  the  PLE  converges  to 

t  i 

exp(-/[F(u)]‘AdF(u,  1)}  for  t  €  [0,  ajCJ5))  Then  we  use  this  result  to  obtain  the 
0 

strong  consistency  of  the  Pl£  under  various  conditions. 

Let  5^(t)  denote  the  PLE.  Then  by  the  continuity  of  G,  f^(t)  is  given  by 


V*> 


1  ,  t(n)  *  0  or  t  e  (-«>,  0), 

^CnFn«n:j)J[nr„»mjJ  *  lj_1- 

*  e  CZn:q-l*  q  *  lf  ***»  t(n)) 

r(n)  , 

£  [nFn(Zn:j^tnFnCZn:j)  ♦  1J*  x(n)  *I,t«  CZ*^,  -). 


Note  that  originally  Kaplan  and  Meier  (1958)  left  the  PLE  undetermined  on  the  sef  • 
(max{Tq,  q  *  1,  ...,  n>,  «). 

We  prove  now  that  ^(t)  converges. 

Theorem  5.1.  Assume  that  Conditions  (3.1),  (3.2)  hold.  Then  there  is  a  set 
e  B,  P{01>  =  1,  such  that  for  all  u>  c 

t  , 

limff  (t)  -  exp(-J[F(u)]  AdF(u,  1)}  for  t  e  [0,  a(F)). 
n+~  0 

t  m  i 

Proof.  Note  that  the  continuity  of  the  function  exp{-/[F(u)  ]  dF(u,  1))  in 

0 

CO,  o(F))  follows  from  the  continuity  of  F(t,  1)  which,  in  turn,  follows  from  the 
continuity  of  G.  Thus  to  prove  the  desired  result  it  suffices  to  show  by  a  standard 
argument  [Chung  (1974),  pp.  132-133]  that: 

t  , 

(5.1)  limffft)  -  exp{ - J[F(u) ] '  aF (u ,  1)},  w.p.l.  for  t  e  [0,  a(F)). 
n*«  71  0 


To  prove  Statement  (5.1),  it  suffices  to  show  by  Lemma  4.4  that: 
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t  , 

(5.2)  lin|  SnFT  ft)  +  f'F(u)2'AdF„(u  1  )|  -  0  -rn  T.  for  t  e  r1,  a(F)). 

n*»  “  0  n  " 

We  prove  now  Statement  (5.2).  Note  that: 

(5.3)  x"1  2  -tnx(l  ♦  x)”1  i  (x  ♦  l)'1  for  x  c  (C,  •). 

By  the  definition  of  fl^t)  and  Inequality  (5.3),  for  t  e  [0,  a(F)): 

a(n,t)  ,  t  . 

-tnfl^t)  i  ^  CnFn^Zn:q)J  *  "  /cFn(u) J ’AdFnG».  "-P-1* 

Thus 

t  , 

(5.4)  lunCirf^t)  +  /[Fn(u)J  dFn(u,  1)  j  2  0,  w.p.l. 

Further,  by  the  definition  of  Mn(t)  and  Inequality  (5.3),  for  t  e  [0,  a(F)) 
a(n,t)  , 

2  J.  [nFn(Zn:q>  *  “  * 

a(n,t)  i  t  i 

£  *  frV*’  *n<*>  « 

-CaFn(t)j_1,  w.p.l. 

Thus,  by  Condition  (3.1): 

t  , 

(5.5)  limCtnFT  (t)  +  /  F  (u)J  <iF  (u,  1)J  s  0,  w.p.l. 

ih-  n  0  n  n 

Statement  (5.2)  follows  now  by  Inequalities  (5.4)  and  (S.5).  || 

Note  that  for  t  e  CO,  o(F)),  lim  max(T  ,  q  ■  1,  ...,  n)  >  t.  Thus 

n*<»  "  _ 

Statement  5.2  holds  regardless  how  we  define  M^t)  on  the  set:  (max(T  ,  q  *  1, 


1/ 


It  follows  from  Theorem  5.1  that  if  we  replace  in  Theorems  3.3,  3.5,  and 
Corollary  3.4,  En(t)  by  f^(t)  and  by  a(F),  the  results  remain  valid.  Thus  we 
obtain  the  strong  consistency  of  the  PLE  under  a  variety  of  Conditions.  In 
particular,  we  obtain  the  strong  consistency  of  the  PLE  when  X^,  are  independent 
r.v.'s,  q  *  1,  2,  ...,  and  under  most  of  the  "traditional"  assumptions. 

Finally  we  note  that  these  results  extend  those  obtained  by  Peterson  (1977) , 
and  by  Langberg,  Proschan,  and  Quins i  (lS.’O) . 


6.  A  Comparison  of  the  Piecewise  Exponential  Estimator  and  the  Product-Limit 
Estimator. 

In  this  final  section  we  point  out  some  differences  and  similarities  between  the 
PEXE  and  the  PLE. 

The  most  obvious  difference  between  the  two  estimators  is  that  the  PEXE  is 
continuous  and  strictly  decreasing  on  [0,  is  t  (n) 5  while  the  PLE  is  a  step  function 
with  jumps  at  the  observed  failures .  Since  in  most  life  testing  situations  the 
survival  function  is  anticipated  to  be  decreasing  smoothly  over  time,  the  PEXE  seems 
the  more  appropriate  estimator  of  a  life  distribution. 

Another  difference  which  favors  the  PEXE  is  its  dependence  on  the  actual 
withdrawal  times  in  each  interval  between  consecutive  observed  failures  (through  the 
total  time  on  test)  compared  with  the  PLE's  dependence  on  only  the  nunber  of  with¬ 
drawals  in  each  of  the  intervals.  The  PEXE  uses  more  information  from  the 
incomplete  data  than  does  the  PLE. 

It  is  clear  from  Sections  4  and  5  that  the  PEXE  and  the  PLE  have  the  same  strong 
(w.p.l.)  limiting,  function.  Also  the  PEXE  has  the  same  weak  limiting  process  as  that 
given  for  the  PLE  by  Breslow  and  Crowley  (1974) .  [A  derivation  of  this  result  is 
forthcoming . 3  Consequently,  finite  sample  comparisons  of  the  PEXE  and  the  PLE  will 
be  important  in  determining  whether  the  differences  cited  in  the  previous  paragraph 
(which  disappear  in  the  limit)  result  in  quantitative  advantages  for  the  PEXE  over 
the  PLE.  Chen,  Hollander,  and  Langberg  (1980)  are  conducting  such  a  study.  They 
assume  that  the  restrictive  assumptions  discussed  in  Section  1  hold  and  that,  in 
addition,  P{X^  >  t)  *  [P{L^  >  t}Jp  for  t  e  [0,  «)>  where  p  is  a  positive  real  number. 
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